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We analyze the Kitaev model on the triangle-honeycomb lattice whose ground state has recently 
been shown to be a chiral spin liquid. We consider two perturbative expansions : the isolated- 
dimer limit containing Abelian anyons and the isolated-triangle limit. In the former case, we derive 
the low-energy effective theory and discuss the role played by multi-plaquette interactions. In this 
phase, we also compute the spin-spin correlation functions for any vortex configuration. In the 
isolated-triangle limit, we show that the effective theory is, at lowest nontrivial order, the Kitaev 
honeycomb model at the isotropic point. We also compute the next-order correction which opens a 
gap and yields non-Abelian anyons. 
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I. INTRODUCTION 

In two dimensions, particles may obey nontrivial braid- 
ing statistics^ However, a direct observation of these 
so-called anyons remains one of the most challenging top- 
ics in physics. Several good candidates have emerged in 
the last years among which the fractional quantum Hall 
effect ' but the braiding of Laughlin quasi-particles has 
still not been performed despite recent proposals based 
on Mach-Zehnder interferometer^. 

Interestingly, at a theoretical level, such exotic excita- 
tions may also arise in spin systems '-''. A simple example 
is provided by the so-called toric code whose elementary 
excitations are known to behave as semions''. Neverthe- 
less, this model is difficult to implement because it is 
based on four-spin interactions which are not easily repro- 
duced in experimental set-ups. A better candidate is un- 
doubtedly the Kitaev honeycomb model' which involves 
only two-spin interactions. Indeed, such a system may 
be realized experimentally in optical lattices either with 
cold atoms''*''''^" or with polar molecules^^. Furthermore, 
in a suitable parameter range, a perturbative low-energy 
effective model of the honeycomb model is the toric code'' 
extended with multi-anyon interactions'". One must 
however keep in mind that, in the honeycomb model, one 
also has fermionic excitations which have to be taken into 
account when braiding anyons' '. The honeycomb model 
has attracted much attention recently^'^'^"''^''''^'' because it 
can additionally be solved exactly via different fermion- 
ization methods (Majorana fermions' or Jordan- Wigner 
transformations^ ' ) . 

One of the most interesting extensions of this model 
suggested in Kitaev's seminal paper has been proposed 
by Yao and Kivelson ^ who considered the same kind of 
model but on the triangle-honeycomb lattice. Indeed, 
in the presence of odd cycles, the system spontaneously 
breaks the time-reversal symmetry and has two topolog- 
ically distinct gapped phases characterized by Abelian 



and non-Abelian excitations. In their study, Yao and 
Kivelson showed that the ground state is a chiral spin liq- 
uid associated to an odd Chern number (note that such 
an exotic state of matter has also been found in another 
spin model^''). Their whole analysis relies on an exact 
treatment of the vortex-free sector (see below for details) 
which allows them to compute the fermionic gap. How- 
ever, as in the honeycomb model, although the ground 
state belongs to this subspace, the low-energy states are 
known to lie in other vortex sectors for a wide range of 
parameters in the Abelian phase. In contrast, near the 
transition point, the fermionic gap in the vortex- free sec- 
tor is smaller than the vortex gap'^'^. The aim of this 
paper is to analyze this low-energy spectrum following 
and extending the procedure developed in Rcf. i j for the 
honeycomb model. 

This paper is organized as follows. In Sec. II, we in- 
troduce the Kitaev model on the triangle-honeycomb lat- 
tice and discuss its symmetries. Section III is devoted to 
the perturbative treatment in the isolated-dimer limit. 
There, we first map the spin model onto an effective spin- 
boson system which is well suited to our analysis. We 
show that the low-energy effective Hamiltonian is related 
to the toric code on the honeycomb lattice although, at 
lowest nontrivial order (six), one only has magnetic-like 
operators. This straightforwardly implies that the low- 
energy excitations are Abelian anyons with a scmionic 
mutual statistics. We also compute the two-spin correla- 
tion functions for any vortex configuration up to order 6 
and check our results for two simple vortex configurations 
(vortex-free and vortex-full) which allow for nonpcrtur- 
bative calculations. Finally, in Sec. IV, we consider the 
isolated-triangle limit ; wc show that the effective low- 
energy Hamiltonian is, at lowest order (one), exactly the 
Kitaev honeycomb model at the isotropic point. The 
next-order correction involves three-spin interactions as 
well as triangular plaquette degrees of freedom. In the 
vortex-free sector, this term is exactly the one studied by 
Kitaev' , which opens a gap, and gives rise to non-Abelian 



excitations. Contrary to the isolated-dimcr limit, one 
cannot diagonalize the effective Hamiltonian for arbitrary 
vortex configurations. Thus, we focus on the vortex-free 
sector and compute the fermionic gap in this hmit, which 
is a check of our perturbative expansion. 



II. MODEL 

We consider the Kitaev model on the triangle- 
honeycomb lattice obtained by replacing each site of the 
honeycomb lattice by a triangle and described by the fol- 
lowing Hamiltonian : 



J — links 



' —links 



(1) 



where a takes values x, y, or z, and links of type x, y, z or 
x', y', z' are illustrated in Fig. 1. In the above formula, 
i and j are the two sites of the a or a' link. Without 
loss of generality', we also consider ferromagnetic inter- 
actions Ja, J'a > 0- This lattice contains six sites per unit 
cell, and two kinds of elementary plaquettes : triangles 
and dodecagons. As in the original Kitaev honeycomb 
model, H commutes with all plaquette operators defined 
as Wp = Yli^p where out(i) denotes the outgoing 

direction at site i with respect to the plaquette p. Note 
that with this definition, plaquette operators have real 
eigenvalues Wp = ±1, whereas, with the convention pro- 
posed in Ref. 7, one has Wp = ±i for odd-loop operators. 

As suggested in Ref. 18, in the following, wc further 
set 

Jx = Jy^J.^ J, and = j; = = J', (2) 

so that the Hamiltonian H respects the symmetries of 
the lattice. It is also time-reversal invariant since it is 
quadratic in the spin operators. However, as explained 
by Kitaev' , the presence of odd cycles (here due to trian- 
gles) breaks this symmetry spontaneously. This symme- 
try breaking may be understood by noting that changing 
the flux of all triangles, i.e., flipping their Wp's, does not 
change the energy so that each eigenstate is, at least, 
two- fold degenerate^''*. 

As in the Kitaev honeycomb model , the Hamiltonian 
H can be mapped onto a free (Majorana) fermion prob- 
lem which allows for an exact solution. Thus, in each 
vortex sector defined by a configuration of the Wp's one 
has a fermionic spectrum. Nevertheless, the low-energy 
states may be given by ground states of other vortex sec- 
tors and, when the corresponding flux configurations arc 
not translation invariant, one is led to solve an impurity- 
like problem. In Ref. 18, Yao and Kivclson numerically 
showed that the ground state of H always lies in the 
vortex- free sector {wp = +1 for all p) and supported this 
analysis by perturbative considerations"". In the follow- 
ing, we shall see that this is verified in the first perturba- 
tive limit we consider (isolated dimers). whereas we did 




FIG. 1: A piece of the triangle-honeycomb lattice which has 
six sites per unit cell. 



not manage to prove it in the other limit. In addition, in 
the isolatcd-dimer limit described in Sec. HI, we compute 
the low-energy spectrum for all vortex configurations and 
give a perturbative expansion of the anyonic gap. 



III. ISOLATED-DIMER LIMIT 

A. Mapping onto an effective spin-boson problem 

Our goal is to perform a perturbative analysis of the 
Abelian phase around the isolated-dimer limit J' 3> J. 
To do so, we shall use the effective spin boson mapping 
introduced in Ref. IL' but, for convenience, let us first 
perform the following rotations : 

(3) 
(4) 

where C is the cycle operator which maps {x,y,z) onto 
(y, z, x). Here, each site is encoded by a cell index i and 
a position index fc = 1 , . . . , 6 inside the cell as shown in 
Fig. 1. Hamiltonian (1) then reads 

i 



a _ 
l.i - 
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^5,i = 


^ "5,i ' 


'^64 = 





With this transformation, the interaction term on each 
dimer (i,j) displayed in cyan in Fig. 1 is simply ai^aj^. 
Denoting 1 1) (||)) as the eigenstate of tr^ with eigenvalue 
of +1 (—1), each cyan dimer can be in four different 
states, 



{|TT),Ui)} with energy -J', 
{|n>,UT>} with energy +J'. 



(6) 



An alternative description of these four states consists 
in interpreting the low-energy (fcrro-magnctic) states as 
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two effective spin states without a quasiparticle and the 
high-energy (antiferro-magnetic) states as two effective 
spin states with one quasiparticle. The energy gap be- 
tween these states is A = 2J' and corresponds to the 
fermionic gap evoked in Sec. II. In the following, we set 
once for all J' = 1/2 or equivalently A = 1. Among the 
possible mappings we choose the following : ' - 

ITT) = l^o),Ui) = |^o), |Ti) = |^i),UT> = l^i), (7) 

where the left (right) spin is the one of the black (white) 
site of the dimer, and double arrows represent the state 
of the effective spin. 

Within this framework, each dimer is reduced to a sin- 
gle site with 4 degrees of freedom [one effective spin 1/2 
and a hardcore boson occupation number (0 or 1)]. Con- 
sidering that each site (• or o, see Fig. 1) of the triangle- 
honeycomb lattice belongs to a cyan dimer, one then has 

^l,=rribl + K) , alo=bl + b,, 

^l-rfibl + K) , al=irt(bl-K), (8) 

where r" 's arc Pauli matrices acting on the effective spin, 
and bl (6^) is the creation (annihilation) operator of a 
hardcore boson at site i which obeys 

[b^,b]]=6,,il-2blk^). (9) 

Hamiltonian (5) then reads 
N 

H = --+Q + n + T+2 + T^2, (10) 
where N is the number of cyan dimers, Q ~ blbi, 

To = -^E0M + ^2;- + 4;i (11) 

i 

t+2=tu^ (12) 

i 

with hopping operators, 

*M = 4/i,i^M^i. (13) 
*2;: = blMr^rli, (14) 
4;' = b[.bs,iry,rl,, (15) 

= -^bl.i+r^2b^■i+r^^l'^l^+r,^. (16) 
*2}+n2 = -^bl-b2A+n2T2,i+ri2' (1'^) 
4f = M,i+r.M^li^ (18) 




FIG. 2: Effective kagome lattice obtained from the triangle- 
lioneycomb lattice by replacing each dimer linking triangles 
by a site. We displayed in gray the "dual" honeycomb lattice 
on which the toric code is defined (see text). 



and pair-creation operators. 



2,i 




(19) 


3,i 


_ u\ V X 
" '^3,i"2,i^2,i^3,i' 


(20) 


l.i 


_ i^t ut y X 


(21) 


2,i+Tl2 




(22) 


3,i 

'2,i+n2 


~ ^"3,i"2,i+Ti2 2,i+Ti2' 


(23) 


1,2-1- Til 

3,i 


^ ^"l,i+Tii"3,i^3,i- 


(24) 



Now, each site is encoded by a cell index i and its position 
inside the cell which takes three values fc = 1,2,3 as 
shown in Fig. 2. Within this formalism, the plaquette 
operators read 

M^A = n^f' VFv = (-l)^'-^'''-n^^ (25) 

iGA iGV 

for triangles and 

Wo = (-l)''^.i''2-+''l.i-"2''l.i-"2+''3--"l''3--"l rf , 

ieO 

(26) 

for the dodecagonal plaquette (which arc hexagonal in 
the effective lattice) located below the cell i (see nota- 
tions in Fig. 2). 

The main interest of this mapping is that the form 
of Hamiltonian (10) is especially adapted to the pcrtur- 
bative treatment developed in the Sec. IIIB. Indeed, 
a key ingredient of our approach which is based on the 
continuous unitary transformations^^ together with the 
particle-number conserving generator^^''^'''^'* is that the 
energy spectrum of the unperturbed Hamiltonian has to 
be equidistant. 
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B. Perturbative analysis of the low-energy sector 

In this section, we used exactly the same method as 
those described in Rcfs. 12 and 25 for the Kitaev honey- 
comb model. Therefore, we skipped all technical details 
and only give here the results of our calculations. 

First, let us note that, in the triangle-honeycomb lat- 
tice, one has 



TV, = 3iVt = 67Vd = 27V, 



(27) 



where N is the number of cyan dimers and s, t, and d 
stand, respectively for sites, triangles (so A^t — + 
A^v)i and dodecagons (or effective hexagons, so iVd = 
A'o). There are thus as many conserved Z2 plaquette 
operators A't -I- as the number N of effective spin 
1/2. This implies that in the low-energy subspace with 
no hardcore boson, the effective Hamiltonian can, in the 
isolated-dimer limit J' 3> J, be written only in terms of 
the plaquette operators, and is thus readily solved. The 
general form of this effective Hamiltonian reads 



H. 



(0) 
off 



Pre I 



(28) 



Pl,---,Pre 



where pi, . . . ,p„ denotes a set of n plaquettes. 

We performed this perturbative expansion of the effec- 
tive Hamiltonian up to order 8 which is the lowest order 
involving multi-plaquette interactions. At this order, the 
constant term is given by 

^ = -l-J^-lj^-lj^ + ^J^. (29) 
iV 2 4 8 64 ^ ' 

The first nontrivial contribution arises at order 6 where 
only hexagonal-plaquctte operators are involved. At or- 
der 8, this latter term is renormalized and triangular- 
plaquette operators come into play. More precisely, one 
has 

Hes =Eo- Ch Wh - ChMM J2 ' (30) 

h h.ti,t2 

where the first sum is performed over all hexagonal pla- 
quettes h (i.e. o) and the second one over triplet plaque- 
ttes made of one hexagon h and two triangles ti and ^2 
(of any kind A or V) adjacent to this hexagon. At order 
8, the coefficients are given by 



h,ti,t2 



33 .8 



(31) 



Thus, at order 6, the spectrum docs not depend on 
the fluxes inside the triangles and this degeneracy is only 
partially lifted at order 8. Note that the signs of Ch and 
Ch,ti,t2 confirm, in this limit, that (one of) the ground 
state lies in the vortex- free sector {wp = +f for all plaque- 
ttes) as conjectured in Ref. 18. In addition, we empha- 
size that triangular-plaquette operators appear by pairs 
which are reminiscent from the time-reversal symmetry 
that the effective Hamiltonian must satisfy (see Sec. II). 



It is interesting to interpret this result in terms of pla- 
quettes and vertex operators. Therefore, one may view 
the effective kagome lattice as a honeycomb lattice where 
each site lies in the middle of the triangles as shown 
in Fig. 2. Within this picture, the plaquette operators 
Wq are interpreted as flux (magnetic) operators, whereas 
VFa,v appear as vertex (electric) operators. In this gauge 
theory language used in the toric-codc model'', our re- 
sults show that, at lowest order, there is no contribution 
of the vertex operators. Thus, the triangle-honeycomb 
lattice, in this isolated-dimer limit does not map onto 
a standard toric-code-like problem. However, the eigen- 
states of the effective Hamiltonian are those of the toric 
code on the hexagonal lattice and, as such, display any- 
onic statistics. To be more precise, one must distinguish 
between electric and magnetic excitations which are lo- 
calized on triangles and hexagons (in the kagome lattice), 
respectively. These two kinds of excitations have mu- 
tual semionic statistics'*' but they individually behave 
as bosons. Finally, let us remark that the gaps of mag- 
netic excitations (order of magnitude J^) and of electric 
excitations (order of magnitude J*) are even smaller than 
the gap in the Kitaev honeycomb model (order of mag- 
nitude J*). This would make an experimental detection 
of anyons in the triangle-honeycomb model even more 
problematic than in the honeycomb model^ 



C. Correlation functions in the low-energy sector 

As in the Kitaev honeycomb model, any correlator in- 
volving an odd number of spin operators vanishes al- 
though the eigenstates break the time-reversal symme- 
try. Indeed, as discussed by Yao and Kivelson^'"", every 
eigenstate is, at least, two-fold degenerate since one can 
flip every triangular plaquette without changing the en- 
ergy but this operation is global. Consequently, as in the 
honeycomb model^ the only nonvanishing correlators 
are products of erf a" on an a dimers. Here, we focus 
on the simplest case involving only one such object, i.e.. 

In the triangle-honeycomb model, one has, a priori, 
nine different functions to consider since the unit cell 
contains nine different dimers. However, with the choice 
of the couplings we made, one only has two different func- 
tions to distinguish : those on "weak" bonds {x, y, z links 
with interaction J) and those on "strong" bonds (x', y', z' 
links with interaction J'). As for the low-energy spec- 
trum, one expects a plaquettc-operator expansion as in 
Eq. (28). We performed the calculation of these two cor- 
relation functions up to order 6 and obtained 



strong 



1-2J^ 



105 



r^Lr^zz_jt(^Wp,+Wp,), 

(32) 

where pi and P2 are the two dodecagonal plaquettes 
shared by the considered strong bond Similarly, 
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since we set J ^ J y ^ wc found for a weak bond (z, j) 



c: 



weak 



8 



63 



J^Wp, (33) 



where p is the dodecagonal plaquctte adjacent to the con- 
sidered bond. 

As can be seen from Eqs. (32) and (33), the presence 
of a vortex is detected at orders 6 and 5, respectively. 
This difference stems from the fact that one analyzes the 
isolatcd-dimer limit for which, at lowest order, (7^^*'°"^ = 

1 whereas Cfr*" = 0. 



D. Checks from Majorana fermions 



1/2 + 

-1/2 



E 



i 



Wt = l 



Wt = -1 



FIG. 3: Spectrum of an isolated triangle with J — The 
eigenstates have quantum number Wt = 1 (left) or = — 1 
(right). 



To check our results we computed exactly the ground- 
state energy in the vortex-free (full) sector for which 
Wp = +1 (-1) for all p using Majorana fermions as de- 
scribed by Kitaev for the honeycomb model' . Following, 
the procedure described in Rcf. "2"), we performed a per- 
turbative expansion of the exact solutions order by order. 
Denoting Cq the ground-state energy per cyan dimer for 



a vortex filling factor ly = 



Number of vortex 
Number of plaquette 



one gets 



-J' 



J2 


J4 


11J6 


U7J^ 


2J' 


32J'3 


128J'5 


^ 8192J'7 




J* 






2J' 


32J'3 + 


64J'5 


8192 



,(34) 
(35) 



Keeping in mind that the number of hexagons is N/3 
and that for each hexagon h, there are 15 triplets 
h, ti,t2, these results are straightforwardly recovered us- 
ing Eqs. (29)-(31). 

One can also check the expression of the correlation 
functions in these vortex configurations. Indeed, the 
Hellmann-Feynman theorem states that 



dJ' 
dJ 



-^Ec:r'=-^r'' (36) 

'^Ecrr' = -2c-r', (37) 



where the sum in Eq. (36) [Eq. (37)] is performed over 
all strong bonds (weak bonds) in the initial lattice. The 
last equalities stem from the fact that for = 0, 1 every 
plaquctte has the same contribution which would not be 
true for other vortex configurations. Using Eqs. (34) and 
(35) and the above relations, one can easily check the 
validity of Eqs. (32) and (33). 



IV. ISOLATED-TRIANGLE LIMIT 

A. Mapping to a spin-boson plaquette problem 

Let us now turn to the isolated-triangle limit J 3> J'. 
These triangles live on the sites of an effective hexagonal 



lattice. For convenience, we use rotated form (5) that was 
already used in the isolated-dimer limit. The spectrum 
of the Hamiltonian of an isolated triangle is made of two 
sets of fourfold-degenerate levels. In each of these sets, 
two levels have eigenvalue Wt = 1 and the other two 
Wt = —1. Setting J = the eigenenergies are ±1/2. 
This information is gathered in Fig. 3. 

As in the isolatcd-dimer limit, we interpret low (high) 
energy states of an isolated triangle as containing zero 
(one) hardcore boson. This hardcore boson degree of 
freedom, which we again denote as b, together with the 
Z2 quantum number wt span a four-dimensional Hilbert 
space. It is then natural to introduce an effective spin 
1/2 to span the full eight-dimensional Hilbert space of a 
triangle. We again denote this effective spin as r (though 
it is not the same as in the other perturbative limit, and 
the same remark holds for b) . The way this effective spin 
is introduced is partly dictated by the operators involved 
in the adjacent dodccagonal-plaquette operators. For the 
triangle (1, 2, 3) of Fig. 1 (and dropping the i index), the 
mapping reads 



-WtSm. (38) 



It is straightforward to check that the above t" operators 
satisfy the usual Pauli matrices algebra. The product 
^2^3 appearing in r^', for example, is the same as the 
one that appears in the dodecagonal-plaquette operator 
having the bond (2, 3) in common with the triangle. It 
is interesting to note that one has to use the operator 
Wt = afa^al to fulfiU the SU(2) algebra. The same 
mapping is used for the other triangles, with 1,2, and 3 
simply replaced by 4, 5, and 6 (see Fig. 1 for notations). 

Since the terms proportional to J in Eq. (5) now read 
—Nt/2 + Q, with Q = J2i where the sum runs over 
the sites i of the effective hexagonal lattice (formed by 
triangles), the last task in rewriting Hamiltonian (5) is 
to find the new form of the terms proportional to J'. A 
simple but lengthy calculation yields 



H = -Nt/2- 



3 



E 

OL—x,y,z Q — links 



E 01 



(39) 



6 



where a indicates a link of type x, y, or z on the honey- 
comb or equivalently brickwall lattice, with the conven- 
tions of Kitaev', as shown in Fig. 4. Furthermore i and 
j are the sites of the effective brickwall lattice on link a, 
and the operators Of reads {Wt.i denotes the triangular- 
plaqucttc operator which is now associated to site i) 



where = = —1, and = 0. 



(40) 
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B. Perturbation analysis of the low-energy sector 

From the above expressions, it is clear that the Hamil- 
tonian can now be written 

H = -^ + Q + To + T+i+T_i+T+2 + T_2, (41) 

where the operators T„ change the number of bosons by 
n and are proportional to J'. With our notations, they 
read 



T 



3 



E E 

a=x,y,z a — links 



(42) 



with 



-'0,i.j 



+2Wt,iWt,j 



(43) 



T-i^tI and T-2^tI 



(44) 
(45) 
(46) 



As in the isolated-dimer limit (where T±i terms are 
absent), with a suitable unitary transformation, such a 
Hamiltonian can be recasted in a unitary equivalent effec- 
tive form HcS which conserves the number of bosons, i.e., 
that commutes with Q. We refer the reader to Ref. 24 
(especially Appendix B), from which it follows that at 
order 2, 



Nt 
2 



- To + [T+i, T_i] + i [T+2, T_2] . (47) 



A tedious calculation then shows that in the low-energy 
subspace with no boson, the effective Hamiltonian (still 
at order 2) has the form 



<^ — f (i + 2.-)-tE E 



J' 



T- T- 



ct Q — links 



+- 



2V3J'' 



(48) 



FIG. 4: A piece of the effective brickwall lattice spanned by 
triangles, together with the notations for x, y, and z links 
(left). Filled (empty) dots represent the triangles made of 
filled (empty) dots in Fig. 1. On the right, we have repre- 
sented one of the second-order three-spin terms that appear 
in Eq. (48). 



The low-energy effective Hamiltonian at order 1 is noth- 
ing but that of the Kitaev honeycomb model at the 
isotropic point, though one should remember that site 
i now also has the Z2 supplementary degree of freedom 
WtA- In Eq- (48), the last (primed) sum is to be taken 
over all possible combinations of three sites i, j, and k 
such that i and k arc nearest neighbors of j, and the spin 
"directions" a, /3, and 7 are such that is an a link, 

{j, fc) is a 7 link, and /3 is the outgoing direction at site 
j of the path {i,j,k) (note that a, f3 and 7 are all dis- 
tinct). For clarity, one such term is illustrated in Fig. 4 
(right). These terms, apart from the plaquette opera- 
tor Wtj, are exactly the ones that arise when switching 
a magnetic field on, in the gapless phase of the Kitaev 
honeycomb model. They open a gap and give proper non- 
Abclian anyonic statistics to the vortices, as detailed in 
Rcf. 7. Actually, one does not know how to diagonal- 
ize analytically H^^ for arbitrary vortex configurations 
even at order 1' and, in particular, how to obtain the 
ground state of each sector. Therefore, contrary to the 
isolated-dimcr limit, one cannot compute the correlation 
functions in this limit. 

Let us remark that the plaquette operators on an ele- 
mentary brick (or hexagon) h, namely Wh = Yii^h r™*''*'' 
are the product of the plaquette operators of the corre- 
sponding dodecagon on the original lattice and of its ad- 
jacent triangles, as follows from Eq. (38). From Eq. (48) 
and the previous remarks, it follows that the vortex-free 
sector contains non-Abelian anyons, which is consistent 
with the findings of Ref. 18. One can also use Kitaev's 
result (see Sec. 8 of Rcf. 7) to obtain the gap in this 
sector, 



A = 6V3 : 



2V3J'' 



= 4J'' 



(49) 



in units where J = , and thus also A = ^^'j if J is 
chosen freely. This value of the gap is consistent with the 
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numerical results obtained by Yao and Kivelson^'*, and 
we also checked it using an expansion of the exact result 
from Majorana fermionization. We also checked that the 
ground-state energies in the vortex-free sector obtained, 
thanks to a Majorana fermionization of Hamiltonian (48) 
or directly of initial Hamiltonian (1), match at order 2. 

Finally, let us mention that at order 1, the ground 
states of Eq. (48) arc such that Wh = 1 for all h. There 
are many such states. The vortex-free state is such a 
state, but any configuration where the six triangles sur- 
rounding one dodecagon are flipped to Wt = — 1 is also 
such a state, since every dodecagon is then surrounded 
by an even number of flipped triangles. From form (48) 
at order 2, we do not know how to prove that the ground 
state is the vortcx-frcc state (we do not even know if this 
is true or if one has to go to higher orders in perturbation 
to prove it). We leave this as an open question. 

V. CONCLUSION 

In this work, we have studied pcrturbativcly the Ki- 
taev model on the triangle-honeycomb model. This 



has allowed us to show that in the isolated-dimer limit, 
the model has low-energy Abelian anyonic excitations, 
whereas in the isolated-triangle limit, the anyons become 
non-Abelian. This picture is consistent with the values 
of the Chern number in each of these phases'^. In the 
isolated-dimer limit, we have furthermore computed the 
low-energy spectrum, as well as the spin-spin correlation 
functions, which both display a plaquette expansion. We 
emphasize that such a computation is not an easy task 
within the Majorana or Jordan- Wigner formalism, which 
are only well suited to study analytically configurations 
of vortices which display translational invariance. 
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